Generalized coarse graining applied to a 4 -theory: A model 
reduction-renormalization group synthesis 
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We develop an algorithmic, system-specific renormalization group (RG) procedure by implement- 
ing a systematic coarse graining procedure that is adapted from model reduction techniques from 
engineering control theory. The resulting "generalized" RG is a consistent generalization of the 
Wilsonian RG. We apply the generalized RG to a 4 field theory. By considering nonequilibrium in 
addition to equilibrium observables, we find that naive power counting breaks down. Additionally, 
a large class of short-wavelength perturbations can drive the system away from both the Gaussian 
and Wilson-Fisher fixed points. 

PACS numbers: 05.10.Cc,05.40.Ca,64.60.Cn,64.60.Ht 



Although it is already known that the renormaliza- 
tion group (RG) is not a black box routine, the purpose 
of this letter is to make it more algorithmic. It is easy 
to be misled into thinking that the RG already is al- 
gorithmic because its key ingredients are coarse grain- 
ing and rescaling the system variables 0, Q- However, 
fully algorithmic implementations of the RG fail for large 
classes of problems because it is not possible to ignore 
the physics of a system and expect to obtain meaning- 
ful results. Capturing the essential physics requires iso- 
lating the appropriate models and the structure of per- 
turbations and uncertainties. Consequently, this process 
is system-specific. Additionally, the scale on which the 
physics is observed must be specified. For instance, for 
bosonic theories the long-wavelength physics is investi- 
gated, while for fermionic system it is the physics near 
the Fermi surface. These considerations suggest that the 
primary obstacles to automation are the model identifi- 
cation and coarse graining processes. In this article we 
present a RG procedure that accounts for these system- 
specific obstacles and apply it to investigate the equilib- 
rium properties of a ^> 4 -theory. 

The perspective adopted in this letter is that the ob- 
stacles mentioned above arise from the identification of 
physical observables. We use control theoretic techniques 
for open systems to systematically identify observables. 
Specifically, consider an open system of the form 



x = f(x) + u, 



(1) 



where x and u are vectors in possibly infinite dimen- 
sional spaces. From a field theoretic or statistical me- 
chanical point of view, without u, equation repre- 
sents the "classical" equations of motion of the system, 
u represents generic driving as well as a possible noise 
to which the system is exposed. If u is generated by 
a continuous stochastic process, this imposes a specific 
structure on the noise. Theories with this restriction on 
u may be mapped to a field theory via Martin-Siggia- 
Rose (MSR)/closed-time-path (CTP) methods UGl In 
these instances, coarse graining is typically dictated by 
conserved quantities and thermodynamic considerations. 
Consequently such systems are locally coarse grained. 



Suppose u is an arbitrary input (e.g. u G L^) to the sys- 
tem. We consider the states or regions in phase space 
that are most accessible via driving to be responsible 
for describing the the essential characteristics of the sys- 
tem. This is analogous to the energy landscape picture 
in statistical mechanics where fluctuations govern which 
states contribute the most to the statistics of the system. 
We use this control theoretic notion of the importance of 
states to specify how to coarse grain and, consequently, 
to generate RG equations. The key step in generalizing 
the RG lies in ascertaining how to coarse grain. 

While equation addresses the effects of perturba- 
tions to the system, it does not allow for the possibility 
of multiscale or constrained observation. We can rem- 
edy this by considering more general open systems of the 
form 



x = f (x) + Bu, 

y = ^x, 



(2) 



where y reflects that only some subspace or, more gen- 
erally, subset of phase space is directly measurable. The 
operators B and C respectively specify the structure of 
how noise may enter the system and which states can be 
measured. The additional structure provides a practi- 
cal way to model real systems and consider experimental 
constraints. The constraint imposed by only measuring y 
strongly influences the relative importance of the internal 
states x and hence coarse graining. For instance, if we 
measure a projected subspace of x over a finite but short 
time horizon, we can expect that those states responsible 
for the transient dynamics will be the most important. 
In this case, conservation laws may play an insignificant 
role in determining how to coarse grain. 

In Q we proposed to coarse grain linear systems based 
on retaining the states contributing most to the response 
of the system to disturbances. Following common prac- 
tice, we coarse grain equation (J2J) based on its lineariza- 
tion about a particular solution with u = 0. In par- 
ticular, to simplify analysis, we only consider lineariza- 
tions about equilibrium solutions. The linearizations are 
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generically described by 



x = Ax + Bu, 
y = Cx. 



(3) 



Associated with equation ® are invariants known as 
Hankel singular values (HSV's) The HSV are non- 

negative real numbers, o- max > Ck > cr m i n that comprise 
the spectrum of the operator W. For systems considered 
over an finite time horizon, tf, consider the positive op- 
erators X and Y, called gramians, that are determined 
by the equations 



^ = AX + Xti + X(0) = 0, 

dtf 

Ay 

— = A ] Y + YA + &C; Y(0) = 0. 
dt f 



W 2 may be factored as 



W 2 



XY. 



(4) 
(5) 

(G) 



HSV's provide a precise measure of the error incurred 
by approximating the effect u has on y with reduced 
order models. The HSV's may be interpreted as supply- 
ing a measure of the importance of the internal states 
x. If W is invertible, it is always possible to find a co- 
ordinate system, called balanced coordinates, such that 
X = Y = diag(cr max , . . . ,a mln ) 0. When equation © 
is transformed to balanced coordinates, the best reduc- 
tions are those that project out the states corresponding 
to small HSV. In other words, the ordering of the HSV, at 
least locally around an equilibrium configuration in phase 
space, specifies how to coarse grain a system. An in depth 
treatment of this material may be found in || 0| . It is 
also sometimes possible to "balance" the full nonlinear 
system [T^ . 

The RG can easily be adapted for HSV-based coarse 
graining. Operator theoretic approaches to the RG 0, Q 
demonstrate that coarse graining in the Wilsonian RG 




FIG. 1: Schematic of shell in k space that is integrated out 
in the (generalized) Wilsonian RG. 



is equivalent to multiplying operators or states by pro- 
jection operators ^tJ- The essense of this work is to 
use HSV's to identify the projection operator. As be- 
fore, suppose that k is a vector index that orders the 
HSV's a K for equation from largest to smallest. A 
generalized Wilsonian RG procedure is obtained by: 1) 
transforming the system to balanced coordinates (about 
an equilibrium solution), 



<£(M) 



R(k, x)4>(x, t)dx, 



so that the partition function takes the form, 

Z[{0}]=/D0exp(-5( g ,{0})) 
= fD4>3exp(-S (&{$}) 



(7) 



(8) 



where g is the original set of coupling con- 
stants/functions, d is the Jacobian from equation 
0, and g is the resulting transformed set of cou- 
pling constants; 2) integrating out K-shells about er m i n 
analogously to wavevector shells; and 3) rescaling k 
and (j) appropriately. An interesting but technically 
challenging variant of this procedure is to integrate out 
tr re -shells instead of K-shells about cr m j n . The remarkable 
feature of this variant is that cr K does not respect spatial 
dimension. For instance, integrating out a single <r K -shell 
may entail integrating out an entire subspace in x-space. 
This seems to provide a natural way to understand how 
three dimensional systems may have regions (e.g. affine 
subspaces) exhibiting one or two dimensional critical 
behavior. The technical challenge lies in rescaling a K . It 
is not clear that rescaling a K will recover the full K-space 
thereby generating a meaningful RG. 

Before applying this procedure to the nonlinear wave 
equation, we apply it to some trivial examples to build 
intuition. We first consider the (driven) diffusion equa- 
tion 



d tC j> = DV 2 (j) + ju. 



(9) 



In this example, B = 7, C = 1, and we take tf — * 00. By 
considering a stable system over an infinite time horizon, 
we only need to solve the Lyapunov equations, 



AX + XA^ + BB ] = 0, 
A^Y + YA + &C = 0, 



(10) 
(11) 



instead of equations Q-©. By taking the Fourier trans- 
form of equations (| 1 1)|> — (| 1 1 1) it is easy to derive that W, 
from equation |(JjJ|, in balanced coordinates is given by 



IF, 



bal 



2D\k\ 2 ' 



(12) 



Here the index, k for the HSV's is just |k|. Thus, for 
the diffusion equation, the most important states are 
those that correspond to small wavevector. Thus, local 
coarse graining is appropriate because the smallest ob- 
servable "fluctuations" are due to the short-wavelength 
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physics. The smallest error is incurred by projecting out 
large wavevectors. We will not complete the RG analysis 
here because the standard RG treatments based on local 
coarse graining are applicable. For a field theoretic treat- 
ment it is possible to obtain the coarse grained action 
from the MSR/CTP formalism llj. however alternative 
formulations may be found in [l], llq . 

As a second example, we consider the (driven) linear 
wave equation 



d 2 4> = v 



V 2 cj) + ju, 



(13) 



By the units of u, it represents a true force acting on 
<f>. This, in addition to the fact that only cf> is the "mea- 
surable" quantity, implies that we have isolated our at- 
tention on 0-based observables. This is choice is a very 
statistical equilibrium and thermodynamic one. We have 
completely neglected tt, the field conjugate to <fi, that rep- 
resents the kinetic contributions to the system. When 
posed as a set of first order equations, equation (|T3*)l be- 
comes 







10" 




"0" 




d t 7T 




v 2 V 2 7 




TT 


(14) 
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By smoothing out the time-cutoff at tt with a damped 
exponential in the integral representation of the solution 
of equations l@J and (JSJ) the problem simplifies to solv- 
ing Lyapunov equations. This smoothing process is also 
known as exponential discounting. With the given form 
of B and C in this problem, we find that the matrix of 
HSV's, W, is approximately given by 



bal 



h/l 



4aw|k| 



I 



2X2: 



(15) 



where ^2x2 is the 2x2 matrix identity, a ~ 1 A/, and <S> is 
the dyadic (algebraic tensor) product. As with the diffu- 
sion equation, short-wavelength physics does not signifi- 
cantly contribute to the response, so locally coarse grain- 
ing is appropriate, we will refrain from treating this ex- 
ample in more detail because the standard RG treatment 
of the wave equation and its 4 nonlinear generalization 
may be readily found in field theory textbooks [T^j . 

The main point of the previous examples is to convey 
that the generalized Wilsonian RG consistently repro- 
duces the results of the standard RG for systems that we 
already intuitively know should be locally coarse grained. 
A great strength of the generalized RG is that it per- 
mits us to tackle less intuitive problems. We apply the 
generalized RG to determine the statistical equilibrium 
properties of a nonlinear wave equation with a partic- 
ular nonequilibrium choice of observables. As will be 
seen, a surprising result is that this choice of observables 
forces us to nonlocally coarse grain. The nonlocality of 
the coarse graining has very interesting implications with 
regard to the resulting induced RG flow. In the remain- 
der of the paper, we then consider some of these novel 
implications. 



The (driven) equations of motion that we are consid- 
ering are 



dt4> = tt + a\U\ 
d t n = V 2 </> + A</> 3 

(3 2 TT 



a 2 u 2 



(16) 



where <f> and tt are real- valued fields. The driving now 
includes generalized forces in addition to "true" forces. 
By expanding around equilibrium solutions of V 2 = 
we find that for each real-space position x, 



B = | °!} ° ) and 

a 2 



C = 



Pi 

02 



(17) 
(18) 



This driving allows for more states in (0, 7r)-phase space 
to be accessible compared to the driving in equation l|14f) . 
This, in combination with the form of y, ensures that 
both <j) and 7r-dependent observables are being consid- 
ered. By using exponential discounting, as mentioned 
earlier, we find that the diagonal operator of HSV's is 
given by 



k|)] 1/2 ®h 



|k|) 



(19) 



Wk does not have the HSV's ordered from largest to 
smallest, so it is not expressed in balanced coordinates. It 
is immediately apparent that the HSV's are large for both 
large and small magnitude wavevector. A heuristic expla- 
nation for this strange result is that for large wavevector, 
tt is a pathologically "fast" variable. However, by driv- 
ing tt with u\ over all wavevector, this permits the fast 
resonances to be excited at large wavevector. The patho- 
logical nature of tt as an observable is analogous to the 
pathological nature of considering £ an observable where 
£ satisfies a Langevin equation [l8| . For this reason, tt 
is a nonequilibrium observable. Furthermore, because 
both the small and large wavelength physics contributes 
strongly to the response of the system, local coarse grain- 
ing cannot be correct. The appropriate coarse graining 
is nonlocal. 

In the case where <p and tt are treated on equal foot- 
ing as observables, which in general may not be the case, 
a i = a 2 = a and (3\ = (3 2 = (3. In the remainder, we 
treat this particular case. Furthermore, without loss of 
generality, we set a = = 1. In this case, equation ltTT)|) 
indicates that the |k| = 1 states are the least important. 
Implementing the second step of the procedure for gen- 
eralized RG involves integrating out k-shells away from 
the |k| = 1 surface. Rather than transform the system 
into the balanced K-coordinates, out of convenience, we 
coarse grain the system in k-space. 

The action for this system, without driving, in Fourier 
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space is given by 



5(g,W) = I 2^p/dk|k| 2 ^(k) + 

h I uLi (e} =1 k,) 0(k 1 )0(k 2 )0(k 3 )0(k 4 ), 

(20) 

where D is the spatial dimension of the system since we 
are only considering the statistical equilibrium properties 
of the system. If we wished to do a full nonequilibrium 
treatment, we should coarse grain the MSR/CTP-action 
for the system. In order to coarse grain, we let </> = </>< + 
4>m + 4>> where <p < is only nonzero for |k| < A, <fi m is only 
nonzero for A < |k| < A -1 , and 0> is only nonzero for 
|k| > A -1 , where A < 1. With this decomposition, the 
path integral measure factors as Dcf) = r D(j> < r D(j) m r D<j> > . 
The RG equations are then induced by integrating out 
4> m and then rescaling the wavevectors and fields. For 
this problem, the rescaling procedure requires that 



*<(k) = 
0>(k) 



^^(A-ik), 
= Z > ^ 2 {Ak), 



(21) 



and p = A _1 k for |k| < A and p = Ak for |k| > A -1 . 
Naive power counting breaks down as a direct result of 
rescaling in the two disjoint wavevector regimes. 

Although we start with a theory where g = (1, A, 0, . . .) 
we can expect that the RG transformations may gener- 
ate new nonlinear terms and that the coupling constants 
may become coupling functions. In fact, g flows towards 
having an infinite number of nontrivial components. In 
particular, the coupling constant A becomes a coupling 
function, A(pi, P2, P3, P4), that may be decomposed into 
the coupling functions {A(j 4 _j)} 4 =1 . Here the notation 
indicates that A( i; j) has i wavevectors with |p„| < 1 and 
j wavevectors with |p„| > 1. If we let A = e~ dl , then to 
first loop order the RG equations for Xu j) (pi, P2, P3, P4) 
are 



di\ 



(v)„ 



+4 1 \ 2 {i J -)J(Pl,P2,P3,P4j 



0(A 3 ), 



(22) 



where 00.4 = D — 4, 014.0 = 4 — 13, 02.2 = —D, 013.1 = 
-2(D - 1), and a 1<3 = -2. Also, J(p%, p 2 , P3, P4) = 
K(Pi + P2) + AT(pi + Pa) + • •• + A'(p 3 + p 4 ), where 
K (q) comes from the one-loop contribution to the 4-point 
vertex and is given by 



2 £ 



dQ 



U<|k|<A-i (2tt) d Q t (Q+W 



(27r)°(l+|q|) 



1 R I D-l D-l \ 



(23) 



x 2 F 1 (D-l 



D-l . 



4|q| 



where B(x,y) is the Euler beta function and 2-F1 is the 
hypergeometric function. Although we do not start with 
a "mass" term in the action (i.e. m 2 </> 2 ), such a term is 
generated by the RG flow. If we denote the mass terms 
for |q| < 1 and |q| > 1 respectively by m?,(q) and m?,(q), 



then the associated RG equations for these terms are 
given by 



dim 2 



Jim < — ^,,o < 



2mj 



SdX 
2(2ir) D 

S D X 
2(2tt) 



D 



0(A 2 ) 
0(A 2 ) 



(24) 
(25) 



The first thing to notice in equation l|22l) is that 
the contribution from tree level, the linear term, indi- 
cates that the coupling functions involving a mixing of 
wavevectors (i.e. i, j ^ 0) are irrelevant. This is the first 
piece of evidence that the theory decouples at small and 
large wavelength. The second piece of evidence for this is 
that K(q) in equation l|23|l diverges as |q| — > 1. Due to 
AT(q)-dependence A acquires through J(pi, P2, P3, P4), 
this divergence seems to act like a barrier in the effec- 
tive coarse grained theory to prevent the large and small 
wavevector physics from mixing. As alluded to earlier, 
this decoupling of the small and large wavelength physics 
is a manifestation of independence of <fi and it as physi- 
cal observables. Rather than being general, we consider 
this decoupling to be special because we are only con- 
sidering a real, scalar c/> 4 -theory. Without driving, this 
model lacks continuous symmetries to be broken. In a 
more general model, the existence of such symmetries 
and the associated set of gauge transformations would 
provide means of coupling <j> and tt. 

It is not possible to ignore the wavevector dependence 
that A acquires because they are relevant for |pj| > 1. 
Specifically, higher derivative perturbations, p n (j) 2 ,n > 2 
and p n (j) m , n > 0, m > 4, in addition to higher order non- 
linearities, cf) n ,n > 4, become relevant when |p.;| > 1. 
However, the couplings at small wavevector, |pi| > 1, 
obey the standard RG equations obtained by local coarse 
graining. Physically this means that just as long as the 
system is not exposed to short-wavelength perturbations, 
the long- wavelength, "thermodynamic" physics will re- 
main robustly observable. If the system is perturbed to 
its large-wavevector regime, then it will flow to a short- 
wavelength fixed point instead of the more familiar Gaus- 
sian and Wilson-Fisher fixed points. This reflects that 
the dynamically faster short-wavelength perturbations 
are able to excite the conjugate field 7r, thereby driving 
the system away from its standard statistical equilibrium. 
Were the conjugate field not accessible to the "noise", 
oi = 0, or not an observable, P2 — 0, this phenomena 
would not have occurred. 

In this letter we have presented a new RG procedure 
and have applied it to a (f> A toy model. This proce- 
dure provides, to the author's knowledge, the first sys- 
tematic means to identify the RG projection operator. 
When both equilibrium and nonequilibrium observables 
are chosen, this RG procedure predicts that naive power 
counting breaks down and that terms that are ordinar- 
ily irrelevant become relevant at large wavevector. Ta- 
ble [I] summarizes these results for the <^> 4 -theory. The 
generalized Wilsonian RG developed here is applicable 
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coarse 
graining 


H < i 

p"0 2 , p n $ m ,n>0, $ n , 
n > 2 m > 4 n > 4 


\p\ > 1 

p (p , p <p ,71 > 0, 0, 

n > 2 m > 4 n>4 


local 
nonlocal 


no no no 
no no no 


no no no 

yes yes yes 



TABLE I: Relevance of perturbations 



to nonequilibrium and heterogeneous systems, finite or 
infinite dimensional systems, and systems with various 
perturbations and uncertainties. Although the RG is 
still formally an uncontrolled approximation, the coarse 



graining is chosen such that the effective, coarsened sys- 
tem is close to the original one. Despite the versatility of 
this method, it is often difficult to analytically determine 
the balancing transformations. However, since there are 
very efficient numerical algorithms for finding balanced 
coordinates, this generalized RG remains a numerically 
useful and practical algorithm. 
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